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We investigate the effects of strong magnetic fields on the QCD phase structure at vanishing 
density by solving the gluon and quark gap equations, and by studying the dynamics of the quark 
scattering with the four-fermi coupling. The chiral crossover temperature as well as the chiral 
condensate are computed. For asymptotically large magnetic fields we find magnetic catalysis, while 
we find inverse magnetic catalysis for intermediate magnetic fields. Moreover, for large magnetic 
fields the chiral phase transition for massless quarks turns into a crossover. 

The underlying mechanisms are then investigated analytically within a few simplifications of the 
full numerical analysis. We find that a combination of gluon screening effects and the weakening of 
the strong coupling is responsible for the phenomenon of inverse catalysis. In turn, the magnetic 
catalysis at large magnetic field is already indicated by simple arguments based on dimensionality. 
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I. INTRODUCTION 

In recent years there has been a growing interest in the 
QCD phase structure in the presence of strong magnetic 
fields, see e.g. m- Such fields may play an important 
role for the physics of the early universe, in compact stars, 
and in non-central heavy ion collisions gl[71 GEES]. 

Despite the rich phenomenology, theoretical predic¬ 
tions are challenging. Starting from QED, e.g. |111414] the 
influence of magnetic fields onto QCD was investigated 
in both model calculations, e.g. |15H29| , such as quark- 
meson, Nambu-Jona-Lasinio models and AdS/QCD, e.g. 
|30H36| . with functional renormalisation group methods, 
e.g. mm, Dyson-Schwinger equations, e.g. jS] IT 1 HI6j 
and in lattice calculations, e.g. BZH53I- 

The importance of magnetic fields for chiral symmetry 
breaking has been pointed out in EH- It has been ar¬ 
gued that chiral symmetry breaking is enhanced due to 
an effective dimensional reduction, the magnetic cataly¬ 
sis. This effect has been linked to an increase of the chi¬ 
ral condensate as well as that of the critical temperature 
T c in model studies. Recent lattice results, g2Ml(5a, 
have shown that while the chiral condensate indeed is 
increased, the critical temperature is decreasing with an 
increasing magnetic field, at least for small enough mag¬ 
netic field strength. This effect has been called inverse 
magnetic catalysis or magnetic inhibition , [S3] . 

Continuum studies have mainly been performed in low 
energy fermionic models, such as the (Polyakov loop en¬ 
hanced) quark-meson- or NJL-model. Hence the reason 
for the discrepancy has to relate to the full dynamics of 
QCD, and in particular the back-reaction of the matter 
sector to the gluonic fluctuations. There have been a 
number of improvements to these model studies to in¬ 
clude QCD dynamics [23] 125H25I [551 [55] . Input parame¬ 
ters of low energy effective models, such as the four-fermi 
coupling, should be determined from the QCD dynam¬ 
ics at larger scales. At these scales they are sensitive 
to sufficiently large external parameters such as temper¬ 


ature, density, or magnetic fields. This has been em¬ 
phasized and used in functional renormalisation group 
(FRG) studies, see [5711501 . The dependence of the four- 
fermi coupling on temperature and magnetic field effects 
including gluon screening has been investigated in the 
recent FRG-work [33) of QCD in strong magnetic fields, 
where inverse magnetic catalysis at small magnetic fields 
and a delayed magnetic catalysis at large fields was found, 
see also [36] for an AdS/QCD computation. 

In the present work we investigate (inverse) magnetic 
catalysis by solving the coupled quark and gluon gap 
equations within the Dyson-Schwinger (DSE) approach 
to QCD, and within a FRG study of the four-fermi cou¬ 
pling based on QCD flows and low energy effective mod¬ 
els. We find magnetic catalysis at large magnetic fields, 
while inverse magnetic catalysis takes place at small mag¬ 
netic fields. 

The present work is organized as follows: The gap 
equations for quark and gluon propagators at finite tem¬ 
perature and magnetic field in two flavor QCD are dis¬ 
cussed in Section [IT] We discuss the dependence of the 
chiral transition temperature T c on the magnetic field as 
well as the magnetic field dependence of the chiral con¬ 
densate. In Section |TTT| the mechanisms behind the phe¬ 
nomena of magnetic and inverse magnetic catalysis are 
evaluated within analytically accessible approximations 
to the gap equations as well as to the dynamics of the 
four-fermi coupling. In this set-up we are also able to 
reproduce the lattice results at eB < 1 GeV 2 . In sum¬ 
mary this provides a complete picture of chiral symmetry 
breaking in the presence of magnetic fields in QCD. 


II. CHIRAL SYMMETRY BREAKING IN 
LARGE MAGNETIC FIELDS 

We investigate chiral symmetry breaking in the pres¬ 
ence of large magnetic fields within a functional contin¬ 
uum approach. To this end we calculate the chiral con- 
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FIG. 1: Quark Dyson-Schwinger equation. Lines with 
blobs stand for fully dressed propagators, vertices with 
large blobs stand for fully dressed vertices. Lines with¬ 
out blobs stand for classical propagators, vertices with 
small blobs stand for classical vertices. 


densate for the two lightest quark flavors and obtain the 
critical temperature T c at finite magnetic field. This is 
done by solving the gap equations for the quark and gluon 
propagator in the presence of a magnetic field using the 
Ritus method |61H67j . The computations are performed 
in the Landau gauge. 


A. Quark and gluon gap equations 

The gap equation for the quark propagator, see Fig. [T] 
depends on the gluon propagator and the quark-gluon 
vertex. The former is expanded about the quenched 
propagator. This expansion has been successfully used 
at vanishing temperature, e.g. [BBi. i59| . and at finite tem¬ 
perature in e.g. [ZOMj, the reliability of this expansion 
has been discussed in m- The quark-gluon vertex is es¬ 
timated with the help of Slavnov-Taylor identities (STIs) 
from the quark and gluon propagators. The systematic 
error of the latter estimate gives rise to the dominating 
systematic error, at vanishing temperature this has been 
investigated in ca, a related upgrade of the vertex will 
be used in a subsequent work. 

The inverse quark and gluon propagators, G q (q) and 
Ga(<1) respectively, read in a tensor decomposition at fi¬ 
nite eB and T 

Gq X (q) = Z q (q) (i73<73 + *7o«)Zo + i r y±q±Z± + M) , 

G~A^{q) = (z n Pf + Z ± P r) q 2 + , (1) 

with P^ v = (g^ -Pifpij'/'pjj) and P± = P - P\\, where 
is the transverse projector. The projection operator 
g^ has the property = Py- The Ritus represen¬ 

tation Eq. ([!]) for the quark propagator is equivalent to 
the Schwinger proper time method, see e.g. Eg. In the 
following we will denote Za = Z\\ and concentrate on the 
Landau gauge, £ = 0. The STIs-induced parametrisation 
of the quark-gluon vertex is introduced as 

r^,p)=7% D J 9 E (g,p), (2) 

with Apf<-f (QiP) discussed in Appendix [a[ The quark gap 


m5wm' = 

FIG. 2: Gluon Dyson-Schwinger equation. The gluon 
line with the yellow dot represents the pure glue loops. 


equation can be written in a compact notation as 

G~\p) = G~l(p ) + Cf ^{gr)G q { q )T v { q ,p)G^{q') , 

Q 

(3) 

with q' = q — p and G q $ as the bare propagator. The 
integration ^ stands for an integration over momenta, 
as well as sums over Matsubara frequencies and Landau 
levels. The gluon propagator can be expanded about its 
pure glue part, 

G^>) = G-^>) + nf(p), (4) 

where we have written the fermionic part of the gluon self 
energy explicitly, while the gluon and ghost loop contri¬ 
butions are contained in G g i ue . The corresponding DSE 
for the gluon propagator within this expansion is depicted 
Fig-0 In the following we consider the back-reaction of 
the vacuum polarisation on the pure glue part as small, 
and approximate 

G~^{p)^G^(jp). (5) 

At vanishing temperature this has been shown to hold 
quantitatively for momenta q > 4 GeV, while for smaller 
momenta this approximation still holds qualitatively with 
an error of less than 20%, see Fig. 6 in |74| . Note 
that for momenta q > 4 GeV the dominant effect of the 
unquenching is the modification of the scales (Aym —> 
Aqcd) and the momentum dependence induced by the 
different /3-functions. This is well-captured with the 
above procedure. In turn, at lower momentum scale the 
non-perturbative mass-gap related to confinement comes 
into play. The magnetic field leads to a shift in the mo¬ 
mentum dependence such as that of the running cou¬ 
pling, as well as (additional) mass-gaps in propagators. 
For both asymptotic regimes (eB —> 0 and eB —> oo) 
these effects are well-captured semi-perturbatively and 
we expect that the approximation © holds well. For the 
intermediate regime we rely on the error estimate at zero 
temperature of about 20% deduced from m- 

The fermionic vacuum polarisation part H^(P) reads 

n /"M = \ il ^ J < y91^)G q {q)T''{q 1 p)G q {q '), (6) 

where the trace includes a sum over the quark flavors. 
Details of this expansion can be found in [g. Here we 
proceed in the lowest Landau level approximation, where 
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FIG. 3: Relation of the quark DSE interaction kernel 
to a 1PI skeleton expansion, which in effect induces an 
effective momentum dependent four-fermi vertex. 


we write down the most general tensor decomposition for 
gluon and quark propagators. Projecting onto different 
tensor compositions, we obtain a coupled set of equa¬ 
tions for the dressing functions of the different tensor 
components. In the next section we will comment on the 
relation of the Dyson-Schwinger equations to other func¬ 
tional expansions and discuss the numerical solutions to 
these equations. 


B. Skeleton expansion 


Before proceeding to the numerical analysis, we discuss 
the standard approximation schemes for the quark-gluon 
vertex used in the Dyson-Schwinger framework from a 
more general point of view. This allows us to connect 
the present ansatze to the approximations used in gap 
equations derived within other functional approaches, 
such as functional renormalisation group (FRG) or nPI- 
approaches. 

DSE studies have made extensively use of the specific 
input for the quark-gluon vertex and the YM-gluon prop¬ 
agator in (A2) and (|A4|) and similar truncations with 


great success. Since the quark and gluon self energy dia¬ 
grams, depicted in Fig.[I]and Fig. [2] contain one bare ver¬ 
tex, the correct renormalisation group behavior and mo¬ 
mentum dependence of the equations must be discussed 
carefully. The truncations to the gap equations (3]) and 
Q can actually be very well motivated from a skeleton 
expansion of the 1PI effective action, which would yield 
similar diagrams as in Fig. [T] and Fig. [2] but with both 
vertices dressed. Fig. [3] serves to strengthen this motiva¬ 
tion as it becomes clear that all approximations should 
encode the correct behavior of the four-fermi interaction, 
which is at the heart of chiral symmetry breaking. This 
allows to consistently reshuffle functional dependencies 
in the interaction kernels of the above equations. 

In turn, the FRG-approach (or nPI effective action) 
can be used to systematically derive gap equations in 
terms of full propagators and vertices respectively, see 
e.g. im Here, we simply note that the 1PI effective 
action can be written as 


r[0] = ^Trlnr[</>] + J d t T k [(j)] — terms, (7) 

where <fr encodes all species of fields, the trace in 0 sums 
over momenta, internal indices and all species of fields in¬ 
cluding relative minus signs for fermions (ip and ip' are 


counted separately), and a logarithmic RG-scale t = Ink. 
The RG-scale in ([T]) is an infrared scale. Momenta 
p 2 < k 2 are suppressed in P fc [cp], and r[^>] = T fc=0 [<?!>]- 
The second term on the right hand side of 0 is a RG- 
improvement term which only contains diagrams with 
two loops and more in full propagators and vertices. To 
see this we discuss the gap equation derived from ([T]). It 
follows by taking the second derivative of |7]) w.r.t. to 
the fields. The first term of the right hand side gives 
the diagrams as in Fig. |T] and Fig. [2] with only full ver¬ 
tices (and additional tadpole diagrams). These diagrams 
can be iteratively re-inserted into the RG-improvement 
term, systematically leading to higher loop diagrams in 
full propagators and vertices. Due to its sole dependence 
on dressed correlation functions such a diagrammatics 
naturally encodes the momentum- as well as the RG- 
running on an equal footing. This also facilitates the con¬ 
sistent renormalisation. Note however that it comes at 
the price of an infinite series of loops diagrams which can 
be computed systematically. Here we take the simplest 
non-trivial approximation which boils down to Fig. ^ and 
Fig. [2] with only full vertices. In terms of the original gap 
equation this leads to the relation 


7 dse 

^qAq 



( 8 ) 


where z^ q is the dressing function of the lPI-quark 
gluon vertex. This immediately leads to the standard 
DSE-dressing in (A2). Moreover, in our numerical study 
the vertices are evaluated at their symmetric momentum 
point. 

Note that, while the ansatz for z is indeed consis¬ 
tent when used in the quark and gluon gap equations, it 
cannot be used in functional equations for higher vertices 
such as the four-fermi vertex. It is already clear from the 
discussion above that a consistent evaluation of renormal¬ 
isation group running and momentum dependence must 
be considered separately for each vertex equation. 


C. Results 

We numerically solve the coupled system of quark and 
gluon functional equations in the lowest Landau level ap¬ 
proximation at finite temperature. This approximation 
is valid in the presence of a clear scale hierarchy with 
eB Aqcd- We use an ansatz for similar to that 
used in Dyson-Schwinger studies, e.g. mm, discussed 
in appendix [A] but adapted for temperature and mag¬ 
netic field effects. 

While at large momentum the influence of temperature 
and magnetic fields is very small, at large temperatures 
and magnetic fields the system is effectively dimension- 
ally reduced and hence the momentum dependencies cor¬ 
responding to the absent dimensions vanish. This can be 
accounted for if we replace by 2|eR| once Q\ < 2|eR| 
and Qq by 2irT for Q 0 < 2ttT as the relevant scale in the 
quark gluon vertex, which is consistent with renormali- 
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FIG. 4: Comparison of the chiral condensate (scenario 
1) for up (continuous lines) and down quark (dashed 
lines) at eB = 12 GeV 2 and eB = 24 GeV 2 . 

sation group arguments. Within this parametrisation we 
are still left to decide what exact momentum scale to 
choose, at which the influence of the external scales T 
and eB is small already. We investigate this question in 
detail in section UlI Al 

The gluon propagator deserves some additional atten¬ 
tion. It is decomposed in different polarisation compo¬ 
nents in the presence of an external magnetic field, see 
e.g. Apart from the splitting into longitudinal and 
transverse components with respect to the heat bath, 
there is an additional splitting transverse and longitu¬ 
dinal to the magnetic field. In the lowest Landau level 
approximation only the polarisation subspace projected 
onto by P^ v = (g^ v — ) receives contributions 

from the quark loop in the self energy, see |B]. Note that 
in analogy to temperature effects, also the other gluon 
components must receive contributions from the interac¬ 
tion with the magnetic field, as gluon and ghost loops mix 
different polarisation components. This is an important 
difference between QCD and QED. From dimensionality 
these contributions are linear in eB at least for asymp¬ 
totically large magnetic fields, leaving aside implicit 13- 
dependencies via the vertices. Their full computation is 
beyond the scope of the present work. Here we investi¬ 
gate the following two limiting cases. 

1. Scenario 1 We simply neglect the screening effect 
of the magnetic field onto those polarisation com¬ 
ponents that feel magnetic effects only through the 
Yang-Mills sector in a QED-type approximation. 
This leads to underestimating the effects leading 
to inverse magnetic catalysis and hence an upper 
limit for T c . 

2. Scenario 2 For the large magnetic fields discussed 
here, the gluon and ghost loops contributions to 
the self energy must have a similar dependence on 
eB as the fermionic part. Since this sector does 
not directly contain charged particles, the effect of 
the magnetic field onto the YM-sector is suppressed 
by powers of the involved couplings. Hence, most 



T [MeV] 

FIG. 5: Comparison of the chiral condensate for sce¬ 
nario 2 at eB = 12 GeV 2 and eB = 24 GeV 2 . 


likely the H-dependence is much smaller than that 
from the fermionic sector. As a limiting case we will 
assume the same magnitude of the self energy for all 
gluon components, which is given by the fermionic 
contributions. With that we overestimate the gluon 
screening effect and obtain a lower limit for T c . 

Both scenarios give consistent limiting cases for the trun¬ 
cation used here. 

As an order parameter for chiral symmetry breaking 
we calculate the chiral condensate as a function of tem¬ 
perature and magnetic field in two flavor QCD in the 
limit of vanishing bare quark masses m u ~ irid ~ 0. The 
Ritus method is not reliable for rather small values of 
qjeB, with qf + 2/3 and —1/3 for up and down quark 
respectively. We expect the lowest Landau level approx¬ 
imation to be a good estimate once eB > 4 GeV 2 (see 
[BJ) which is also the regime where the approximation 
(15]) works well for vanishing temperature. 

The numerical computation is very demanding in the 
vicinity of the phase transition due to the diverging cor¬ 
relation length. This translates into a numerical error 
in the critical temperature indicated by the error bars 
in the plots. Fig. [4] and Fig. [5] show the up- and down- 
quark condensate for different values of eB. The inverse 
magnetic catalysis effect described in (Mills] is evident. 
While the chiral condensate still rises with the exter¬ 
nal field in the low temperature limit, the transition be¬ 
tween chiral broken and symmetric phase drops. This 
signals inverse magnetic catalysis as observed on the lat¬ 
tice, |47l {48} . Furthermore the phase transition, which is 
second order at zero magnetic field turns into a crossover 
with growing eB, even for vanishing bare quark masses. 
This can be understood as magnetic screening: the mag¬ 
netic field effectively serves as an infrared cutoff, which 
inhibits an infinite correlation length. 

In the present computation in two-flavor QCD, an even 
more intricate effect is observed. Up and down quarks 
come with different electric charges, therefore the pres¬ 
ence of a strong electromagnetic field breaks isospin ex¬ 
plicitly. This results in a non-degenerate chiral phase 
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FIG. 6: Comparison of the chiral condensate at zero 
bare mass and at a finite bare quark mass of m u = 
rrid = 10 MeV at eB = 4 GeV 2 in scenario 1. 


transition for the two flavors. Because gluons travel 
through a medium filled with both virtual up and down 
quarks, isospin breaking effects the self interactions of 
the quarks, which leads to interference between the chi¬ 
ral transitions of the two flavors as seen in Fig. [4] and 

Fig- 111 

This interference can be interpreted as follows. Virtual 
quark fluctuations contributing to the gluon screening are 
suppressed in the chiral broken phase by the quark mass. 
Since the down quark undergoes the chiral phase transi¬ 
tion already at lower scales, its fluctuations are suddenly 
enhanced due to the vanishing mass in the symmetric 
phase. The up quark, while still in its chirally broken 
phase, is drastically effected by these enhanced fluctua¬ 
tions, which lead to reduction of the up quark condensate 
even below the real phase transition. 

It can be seen from Fig. [4] and Fig. [5] that this effect 
is more prominent in scenario 2, which should come as 
no surprise, as the coupling of the magnetic field to the 
gauge sector is probably overestimated here. Neverthe¬ 
less the isospin induced chiral transition substructure is 
observable in the limiting scenario 1 as well, which is a 
strong indication of its validity. Therefore this important 
physical effect might be observable in lattice calculations, 
as well. In [48l [49] the averaged chiral condensate was 
investigated at finite quark mass. However when we in¬ 
vestigate the chiral transition at a bare quark mass of 
10 MeV we find that the interference effect is completely 
masked by the crossover behavior as can be seen in Fig. [6] 
Note that here the unregularized condensate at finite bare 
mass is plotted, hence the offset between the curves. 

In analogy with lattice calculation we define T c at the 
inflection points of the curves shown. In Fig.[7]and Fig. [8] 
the obtained values for T c for the limiting cases described 
by scenario 1 and 2 are shown. The two curves give lower 
and upper limits for T c , as discussed before. The chiral 
transition temperature is decreasing for a large range in 
eB before it seems to saturate for intermediate values in 
both scenarios. At very large fields it rises again. 



FIG. 7: Critical temperature obtained from scenario 1 
for up quark, down quark and from the flavor averaged 
condensate. 


In accordance with our previous discussions we see that 
the up and down quark chiral transitions do not coin¬ 
cide. The transition temperature from the flavor aver¬ 
aged quark condensate is given in Fig. [7] and Fig. [8] as 
well. As can be seen from Fig. [4] and Fig. [5] the tran¬ 
sition temperature of the flavor averaged condensate is 
essentially determined by the up quark. 

Both scenarios give estimates for the chiral transi¬ 
tion temperature, which differ only quantitatively. Sce¬ 
nario 1, which underestimates the magnetic field effects 
in the gluon sector extrapolates to a critical tempera¬ 
ture at eB = 0 between 170 — 210 MeV with a turning 
point between catalysis and inverse catalysis of about 
eB ss 30 GeV 2 . On the other hand scenario 2 gives T c 
at zero magnetic field of about 140 — 165 MeV with a 
turning point slightly higher than in scenario 1. This is 
in accordance with the fact that scenario 2 overestimates 
the gluonic sector, which is the source of the inverse catal¬ 
ysis effects. At B = 0 the chiral phase transitions for up 
and down quark coincide. While the continuous lines in 
Fig- 0 and Fig. [8] are obtained from a fit with a simple 
quadratic polynomial, reflecting the turnover behavior at 
large fields, these should not be mistaken as extrapola¬ 
tions towards zero. Furthermore the computations have 
been performed in the lowest landau level approximation. 
This leads to an uncertainty of about 10% for B smaller 
than 10 GeV 2 , while the qualitative behavior is not ef¬ 
fected, as discussed in 00. In the following section we 
will see that the behavior of T c at small B is steeper than 
just quadratic. 

It is well known that within approximation schemes 
such as the one discussed here, relative fluctuation scales 
are usually well accounted for, whereas absolute scales 
have to be fixed. The position of T c at eB = 0 gives us 
the possibility of identifying absolute scales and allows to 
adjust our truncation. We will not be concerned about 
matching the exact scale of T c at zero magnetic field with 
the lattice, moreover we will investigate the mechanisms 
behind the B — T phase structure in greater detail. We 
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FIG. 8: Critical temperature obtained from scenario 2. 


will discuss the issue of scales in the following sections. 


but a manifestly renormalisation group invariant approx¬ 
imation of the above Dyson-Schwinger equations, see the 
discussion in Section |IIB| It includes only dressed ver¬ 
tices. In appendix |A| we discuss how the interaction ker¬ 
nels can be related in both pictures. The 1PI quark-gluon 
vertex is parametrized as 

r gA 9 (9 2 ) = z a 2 ^a s (q 2 )^ , (10) 

The gluon propagator is transversal due to the Lan¬ 
dau gauge, and we allow for a gluonic mass via thermal 
and magnetic effects. M(py) is a function that is ap¬ 
proximately constant in the IR but falls of rapidly for 
Py > 2|eR|. Hence, if we are interested in M{ 0) = Mir 
we can write, dividing the equation by its trivial solution, 


1 - 4:TT 2 C F T y 


M k + <l\u 


III. ANALYTIC APPROACHES 

In the present Section we are specifically interested in 
the mechanisms at work in magnetic and inverse mag¬ 
netic catalysis. To that end we discuss approximations 
to the quark gap equation in Section |III A| as well as to 
the dynamics of the four-fermi coupling or quark scat¬ 
tering kernel in Section |IIIB| that allow for an analytic 
approach to chiral symmetry breaking. While the quark 
gap equation can be straightforwardly reduced to an an¬ 
alytic form from that used for the numerical study, the 
four-fermi coupling is studied in a renormalisation group 
approach to QCD, that reduces to an NJL-type model 
for low momentum scales. 


A. Quark gap equation 

The mechanisms behind the phenomena observed in 
our numerical study can be analyzed within approxima¬ 
tions detailed below, that allow for an analytic access. 
These approximations to the gap equation have been in¬ 
troduced in m for QED, and can be extended to QCD at 
finite temperature. The self-consistent Dyson-Schwinger 
equation for the mass functions reads in lowest Landau 
level approximation with zero bare mass 


x 


^ ots exp {—x/2\eB\) 
q^ b + x + U( x ,q lb ) 



%b + X J 


= 0. 

( 11 ) 


In (111 we have introduced q^b = (q3,2mrT) and q\\j = 
(g 3 ,Z7rT(n-|-l/2)). Chiral symmetry breaking is realized 
once a solution M 2 R > 0 exists. Due to the shape of M(q) 
and the exponential factor in (11), the integrand only 
has support for x < 2|e£>|. In the following we carefully 
investigate the ingredients to this self consistent equation 
and the physical mechanisms, which are responsible for 
the intriguing behavior seen in the previous section. 


Due to the finite support of the integrand, the mo¬ 
menta running through the vertices are comparable or 
smaller than the relevant dimensionful quantities eB and 
T 2 . Note that in our numerical study we have used an 
ansatz for the quark gluon vertex, that includes generic 
eB and T dependencies. Here we utilize the fact that 
the running of a s is dominated by the temperature and 
magnetic field scales. We resort to a simple ansatz for 
a s(Q 2 /^Qcn) based on the analytic coupling a SjHQ sug¬ 
gested in (SDJ HU], see [82] for an investigation within the 
present context. This coupling yields a linear potential 
such as seen in the heavy quark limit. 


<x s (z) = a s , HQ (z) r m (z ), (12) 


... , „ ~ vf M (9||) Tr ( A (sgn(eR))7y 1 7||) 

M(P " ) = 4 ’ Cf J-- 

t ( k l \ CAt) 

J °‘ exp (~2^i) WTmp-y (9) 

k± 

Here = Tj2 n J dgj|/(27r) 3 and A(s) = (1 + scr 3 )/ 2. 
The quark gap equation d9| is obtained from a skeleton 
expansion of the effective action, e.g. eg, and is nothing 


where 


u s , uq { z ) = 


1 2 2 - 1 


/3 0 z 2 log(z 2 ) ’ 
with /3 q = (33 — 2Nf)/Y2-K and 


* 2 = 


Afj2 eB + Xt(2ttT) z 
X 2 

7V qcd 


(13) 


(14) 


with coefficients At, A b, which are of order one. These 
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coefficients determine the point at which eB or T domi¬ 
nate momentum scales. For the relevant magnetic fields 
and temperatures the running of the coupling with tem¬ 
perature is very small compared to the running with eB. 
We use an ansatz for the infrared behavior of the vertex, 
which is parametrized in riR. Here we use 


ri R {z 2 ) = 


(z 2 + b 2 ) 2 


1 + 


b 2 


(15) 


which scales with oc z 4 for z —> 0, and approaches unity 


in the perturbative regime. Eq. (121 reproduces the cor¬ 
rect behavior of the full quark gluon vertex in (101. We 


leave b and c as parameters which allow us to model the 
infrared behavior of the quark gluon vertex. Our ansatz 
for (15) is motivated from the quantitative renormalisa¬ 


tion group study of quenched QCD in m, which we use 
to determine b and c. We get 


6=1.50, c = 7.68 , 


(16) 


from the fit to Fig. 4 in E2- 

Furthermore we discuss the gluon self energy in the 
presence of magnetic fields at finite temperature in this 
simplified setup. It is important to notice that we can 
facilitate our calculations by the following argument. The 
function on the right hand side of (11) is a continuous real 
function of Mir and approaches +1 as Mir —>• oo. Hence 
it is sufficient to check whether the expression is negative 
for Mir = 0, because then it had to pass through zero at 
some point, which means that a solution exists. 

The gluon self energy receives two important contri¬ 
butions. The first is through the appearance of fermion 
loops, which are also present in an abelian calculation. 
The fermionic self energy part in lowest Landau level ap¬ 
proximation with Mir = 0 factorizes 

H$r (p) = aeB exp (-p 2 j2eB)H^(p ilf T) . (17) 

Contracting with P^ in the Landau gauge, we can write 
the second term as 


H/(P||,T) = -8tt 2 3-2(1 -p\/p 2 ) 


1 


x dx 


x(x — 1) 


S (?1 + (2 n) 2 (n + 1/2)2 + x{1 _ x)/ t 2 Y 

(18) 

where we defined r 2 = T 2 /p |. The function can be evalu¬ 
ated numerically and is very well described by the simple 
function 


H/(ph j T) = (1/27t) 3 — 2(1 — pyp 2 ) 


1 


1 + (47t 2 /3)t 2 

(19) 


Eq. (171 and Eq. (19) state that the relevant contribu¬ 



FIG. 9: Analytic calculation of the critical temperature 
for the chiral phase transition. The bands indicated 
correspond to At = 1 and A t = 0. Arrows indicate the 
direction from At = 1 to At = 0. 



FIG. 10: Comparison of the critical temperature ob¬ 
tained with our full numerical procedure to the simple 
analytic estimate for Ar = 1.1, At = 1 and n = 1.19. 


tions to the self energy stem from p/ ~ 2 eB and p| a: T 2 . 
Similar as before, the influence of the magnetic field onto 
the Yang-Mills sector is not easily accounted for. Here we 
focus on the abelian-like part of the gluon self energy. As 
we have investigated before numerically, this is qualita¬ 
tively correct and we will use Eq. (14) to account for the 
correct scales. It is well known from Dyson-Schwinger 
studies [83], that approximations similar to this semi- 
bare vertex ansatz underestimate the strength of chiral 
symmetry breaking, due to the negligience of important 
tensor structures in the vertex, especially those struc¬ 
tures that break chiral symmetry explicitly m- In order 
to compensate the overall weakness of the interaction, we 
allow for a phenomenological parameter n in front of the 
integral in Eq. Ca¬ 
using our simple ansatz we can investigate chiral sym¬ 
metry. In Fig. [9] a family of solutions to Eq. (111 is shown 
for various values of Xb and At, using the ansatz de¬ 
scribed above with k = 1.2 for the two upper curves and 
k = 1.4 for the lower curves. The choice of n is for better 
























B. Four-fermi coupling 
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FIG. 11: Diagrams contributing to the renormalisation 
group flow of the four-fermi coupling. 


visualisation only, as the curves can be shifted up and 
down using this parameter. 

The observed behavior agrees with that in our numer¬ 
ical study. It can be seen from Fig. [9] that for small eB 
inverse magnetic catalysis is present, while at large eB 
the the critical temperature rises again with the magnetic 
field, with 


Tc{B/ Aqcd ~t oo) c* y feB , (20) 

as one would anticipate from dimensional considerations. 
This behavior is universal for all A b and A t- We see that 
the choice of A^ effects the position of the turning point 
of the chiral phase boundary. 

With the present analytical considerations the numer¬ 
ical results in Fig. [7] and Fig. [8] are readily explained: 
they roughly correspond to A b ~ 1, which explains the 
relatively large value of eB at the turning point. We see 
that already small changes in As have a huge effect on 
this quantity, see Fig. [9] 

In Fig. [TO] we have plotted the analytic result with 
A b = 1.1, At = 1 and k = 1.19, which agrees well 
with the numerical results from scenario 2. Based on 
the present work we estimate that A b « 2 — 3 is a realis¬ 
tic choice for the B-dependence of the running coupling, 
as in our numerical study quark and gluon propagator 
turn into their corresponding B = O-propagators at this 
momentum scale. 


For a further analytical grip we also resort to a low 
energy effective theory point of view: integrating-out the 
gapped gluons leads to an effective four-fermi theory, that 
is initialized at about the decoupling scale of the glue 
sector of A k 1 GeV. Previously there have been phe¬ 
nomenological approaches in low energy effective models 
to include QCD dynamics as the source of the inverse 
magnetic catalysis effect |23J [25, [55]. From the point of 
view of the FRG for QCD this can be seen as follows 
M EDI EJ US I84H86] : At a large momentum scale k 
QCD is perturbative, and the 1PI effective action in 
([7]) is well-described perturbatively. A four-fermi cou¬ 
pling is generated from the one-loop diagrams (in full 
propagators and vertices) encoded in 0 , the related di¬ 
agrams are depicted in Fig. HH In the present discussion 
we have dropped diagrams that depend on the qq — AA 
vertex, qqqq — AA-vertex and the qqqqqq-ve rtex. Fur¬ 
thermore we assume a classical tensor structure for the 
qAg-vertex with a coupling ^/47ra Sj fc, and only consider 
the scalar-pseudo-scalar four-fermi vertex 

IWfermi [q, q, B] = \ Tg$ m f q$ , (21) 


with the scalar-pseudo-scalar tensor structure 

r k% d im = A k [6 ij 6 l mS ab 6 cd + (rysMi ^)im(r n ) ab (r n ) cd ] • 

( 22 ) 


The four-fermi term in (211 can be viewed as the inter¬ 
action term of a NJL-type model. Within the approxi¬ 
mation to QCD outlined above the flow of the four-fermi 
coupling, d t Xk , has the form 

dtXk = — k 2 — \kOt s ,kF\a a {G qi GA) 


at 


--%TF a *(G q ,G A ), 


(23) 


The present analysis reveals the following mechanism: 
The gauge sector acquires a B-dependence through the 
feedback of the fermionic sector. This dependence is 
responsible for the phenomena called inverse magnetic 
catalysis, as has been also observed recently in a FRG- 
study within QCD, |T3]. This also explains why it cannot 
be seen in model calculation without explicit QCD input. 
From Eq. (Ill, Eq. (121 and Eq. (191 we see that the 


gluon screening and the running of the strong coupling 
(both by thermal and magnetic effects) are competing 
with the generic fermionic enhancement of chiral sym¬ 
metry breaking in a dimensionally reduced system. We 
see from Fig. [9] that at small magnetic field screening ef¬ 
fects dominate the behavior of the fermionic self energy, 
while at asymptotically large fields, thermal fluctuations 
are negligible and hence eB, as the dominating scale, 
drives the phase transition towards higher T c (magnetic 
catalysis). 


with positive coefficients F\, F\ as , F as . The respective 
diagrams are depicted in Fig. ED The different classes 
of diagrams in Fig. 11 depend on combinations of gluon 
and quark propagators, Ga and G q respectively. 

The four-fermi coupling A& in two-flavor QCD at 
T = 0 has been quantitatively computed (including 
its momentum-dependence) in quenched QCD with the 
FRG in [75 , and in a more qualitative approximation 
(without its momentum-dependence) in fully dynami¬ 
cal QCD in [73]. The respective results are depicted in 
Fig. [12] As expected, the couplings have a similar de¬ 
pendence and maximal strength. However, the slope of 
the coupling in the qualitative computation in the peak 
regime relevant for chiral symmetry breaking is bigger for 
the qualitative computation. This can be traced back to 
the missing momentum-dependencies, whose lack artifi¬ 
cially increases the locality in momentum space and in 
the cutoff scale. Hence, guided by the experience gained 








9 



k [GeV] 


FIG. 12: Scalar-pseudo-scalar four-fermi coupling in 
the vacuum, T = 0, B = 0, computed with quantita¬ 
tively reliable QCD-flows in quenched QCD, ca. and 
with qualitative full QCD flows, m- 


in the DSE-computations we expect the slope to play a 
large role and we shall use the quantitative quenched re¬ 
sults for Afc and a s in our present computations. We shall 
further comment on the differences in the next Section. 

For large cutoff scales k the propagators approach the 
classical propagators. The current quark mass at these 
scales is negligible and only the cutoff scale is present, if 
temperature and magnetic field are considered small rel¬ 
ative to the cutoff scale. Then the dimensionless Fs are 
simple combinatorial factors. For optimized regulators, 
m , they are given as 


F\ = 4IV C , 


F Xa , = 12 


N 2 C - 1 
2 N c 


F„ 


3 97V 2 - 24 

16 W c ’ 

(24) 


in the vacuum, see e.g. EMM for more details. For 
small enough cutoff scales k the gluonic diagrams decou¬ 
ple due to the QCD mass gap. In the Landau gauge 
this can be directly seen with the gapping of the gluon 
propagator. For T = 0, B = 0 this entails 


p 2 G A (p 2 < A 2 ) (xp 2 /r 


gap ' 


(25) 


with A « 1 GeV. We emphasize that (25) only reflects 


the mass gap present in the Landau gauge gluon prop¬ 
agator, the gluon propagator is not that of a massive 
particle, see e.g. [88] ■ For momentum scales p 2 < A 2 this 
approximately leaves us with an NJL-type model with 
the action 


Lnjl \Qi Qi B] — / Q “b Lf our _f erm j[ q, cj, .B], ( 26) 


with the scalar-pseudo-scalar four-fermi interaction de¬ 
fined in (211. In the presence of a magnetic field this 
model including fermionic fluctuations has been inves¬ 
tigated in [38, within the FRG. Here we shall use the 
respective results within the lowest Landau level approx¬ 
imation. Then T c shows an exponential dependence on 


the dimensionful parameter eB 


T c = 0.42A exp 


/ 


V 


2t r 2 


n c Xa E Wf eB \ 

f 


(27) 


The well-known exponential dependence of T c on the 
four-fermi coupling Aa already explains the large sensi¬ 
tivity of the scales of magnetic calatysis and inverse mag¬ 
netic catalysis to details of the computation. Eq. ( f27| ) is 
valid for large magnetic field and for A -C m 2 ap , that is 
deep in the decoupling regime of the gluons. An estimate 
that also interpolates to small magnetic fields is given by 


T c = 0.42A exp — 


1 


caAa 


(28) 


with 


N, 


Ck{ 


( B ) = 2 ^ 2 1 +cik2 


with ci = 3 , 
(29) 


where ci has been adjusted to reproduce T C (B = 0) ~ 
158 MeV. While Eq. (291 resembles a lowest Landau level 


approximation, it is actually an expansion in B. Using 
this ansatz we can describe the behavior of the phase 
transition on scales below 1 GeV 2 qualitatively, while the 
B = 0 limit is fixed. 

It is also well-known that for k m gap the flow of 
the four-fermi coupling is driven by the gluonic diagrams 
summed-up in F as : for large scales we can set Afc>m ga . p ~ 
0. The gauge coupling is small, a s ,/c>m gap *C 1 and the 
flow gives A k oc a 2 . This entails that the diagrams with 
four-fermi couplings are suppressed by additional powers 
of a s , and the four-fermi coupling obeys 


d t X 


glue,fe — 


^S,k 

~w 


F aa ( Gq , Ga) 1 


(30) 


where the subscript ’glue’ indicates that the flow is driven 
by glue fluctuations. As discussed before, for k 3> m gap 
we have classical dispersions for quark and gluon, and 
the diagrammatic factor F as is a constant, see (241. The 


strong coupling a s k h as the form (12) with z oc k. Inte¬ 
grating (301 with (121 gives 


\iue,fc oc -^F a ,(G g , Ga) ■ 


(31) 


where an estimate for the H-dependence of the gluonic 
diagram in F as is given in Appendix |B| 

At vanishing magnetic field A gluo ^ agrees well with the 
full result for the four-fermi coupling in [75] for k > 2 
GeV, see Fig. 13 Below k « 2 GeV, A glue ,fc is increas¬ 


ingly smaller than the full scalar-pseudo-scalar four-fermi 
coupling in quenched QCD. In this intermediate range, 
where all diagrams contribute, we write the resulting cou- 
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FIG. 13: Scalar-pseudo-scalar four-fermi coupling at 
T = 0, B = 0 computed with quantitatively reliable 
QCD-flows in quenched QCD, da, in comparison to 
Agi U e computed from (311. 


pling within a resummed form that captures already the 
fermionic diagram proportional to F \, 


Afc = 


A k 


1 0^ Afc 


with Ck = [ dk'k'F\{G q ). (32) 

Jk 


The resummed form in 


already reflects the matter 
part of the flow in (|23| which is the term proportional to 


dtXk- The other terms add up to 


&tXk — (1 C-kXk) ( A k^s,k 



For CfcAfc <C 1 the flow of A k boils down to (30|). 
CfcA —► 1 the flow in (33) tends towards zero. 


(33) 

For 
n this 


regime the four-fermi coupling grows large and the mat¬ 
ter flow dominates. Hence, for the present qu alita tive 
analysis we simply identify A with the glue A glue) (31), up 
to a prefactor, 


A k — Z\\, 


A ''glue, k • 


(34) 


The prefactor Z\ accounts for the fact that we have used 
results of quantitative QCD-flows m for the strong cou¬ 
pling which also includes wave function renormalisations 
for the quarks. In the current model considerations with¬ 
out wave function renormalisation and further simplifica¬ 
tions this has to be accounted for. For the same reason 
the normalisation 0.42 A related to a four-fermi flow with 
an optimised regulator has to be generalised. Moreover, 
the prefactor c\ t k is the integrated four-fermi flow already 
present in (281 up to an overall normalisation accounting 


for the model simplifications. We choose 


Ck{B) = c 3 Ck{B ), and 0.42 A —>• 0.42 A exp (c 2 — c 3 ) , 

(35) 


and arrive at 


T c = 0.42Aexp ( — 


caAa 


c 2 


(36) 


with ca as given in (291 and A in (341 and (311. Note 


that the parameter c 3 has dropped out. Its value can be 
adjusted to achieve a quantitative agreement of (311 with 
the QCD result in [75] with 


C3 = 


1 

2Za 


(37) 


where the factor 1 fZ\ simply removes the mapping factor 
adjusting for the missing wave function renormalisations 
in the model computation. This quantitative agreement 
strongly supports the reliability of the approximate solu¬ 
tion to the flow equation given by (321 in the intermediate 


momentum regime that is of importance for the current 
considerations. The remaining parameters are fixed as 
follows, 


Za = 2.2 , Cl = 3, c 2 = 1.4 . 


(38) 


The parameter c 3 has already been adjusted to meet 
T C (B = 0) « 158 MeV, see (28) and (29). The parameter 
c 2 re-adjusts the overall scale 0.42 A —► 0.42 A exp c 2 = 
1.7A. As already discussed above, it depends on the reg¬ 
ulator and the approximation at hand. It reflects the 
dependence on the renormalisation group scheme. Simi¬ 
larly to ci it is fixed with T C (B = 0) « 158 MeV, and is 
a function of the overall normalisation of the four-fermi 
coupling Z\. The latter is the only free parameter left. 
In ( [38| we use the value that reproduces the lattice re¬ 
sults, see Fig. [14] We emphasise that no other parameter 
is present that allows to shift the minimum in T c , the lat¬ 
ter being a prediction. 

Obviously, the effect seen in our numerical and ana¬ 
lytic DSE-study, is also present in the analytic approach 
to the dynamics of the four-fermi coupling, including a 
direct grip on the underlying mechanisms. We see that 
the non-monotonous behavior, i.e. the delayed magnetic 
catalysis, [33! [52], is already present at smaller scales 
compared to Fig. [7] and Fig. [8] while the lattice results 
are reproduced. 

In turn, for asymptotically large magnetic field, the 
critical temperature runs logarithmically with B , 


T C {B/ Aq CD —> oo) a IiiH/Aqcd , 


(39) 


related to a double-log-dependence on B of the expo¬ 
nent. Due to the qualitative nature of the approximation 
of the H-dependence of the gluon propagator it cannot be 
trusted for asymptotically large B. Indeed, (391 has to be 


compared to (201 within the analytic DSE-approach pre¬ 


dicting a square root dependence. Note that in the latter 
computation the quark vaccum polarisation is included 
selfconsistently at large B even though the backreaction 
on the pure glue loops in Fig. [2] is neglected. Still this 
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FIG. 14: Comparison of the chiral transition temper¬ 
ature obtained within the simple mean field NJL es¬ 
timate Eq. (141 to the lattice results of [47] (see their 
Fig. 10). 


indicates the validity of the square root dependence, even 
though a definite answer to this question requires more 
work. 


C. Discussion of scales & mechanisms 


With the findings of the last two sections we have 
achieved an analytic understanding of the mechanisms 
at work. The decrease of T c for small magnetic fields, 
the increase of T c for larger fields, as well as the related 
magnetic field regimes can now be understood. In partic¬ 
ular this concerns the magnetic field R min , where T c (B min ) 
is at its minimum. This is the turning point between in¬ 
creasing and decreasing T C (B). 

Magnetic catalysis relates to the dimensional reduction 
due to the magnetic field in diagrams with quark corre¬ 
lation functions leading to an increase of the condensate. 
At finite temperature the catalysis due to the dimen¬ 
sional reduction is accompanied by a thermal gapping of 
the quarks that counteracts against the magnetic cataly¬ 
sis effects. In total this leads to a rise of both, the chiral 
condensate and the critical temperature, if the magnetic 
field dependence of the involved couplings is sufficiently 
small. As the magnetic field also sets a momentum scale 
of the physics involved, this scenario holds true for suf¬ 
ficiently large magnetic field strength eR/A 2 CD 1, 
where the R-dependence of the couplings can be com¬ 
puted (semi-)-perturbatively. This explains the regime 
of delayed magnetic catalysis. 

The above discussion of the standard scenario already 
entails that rapidly changing couplings are required for 
a decreasing T c . The couplings involved are the scalar- 
pseudo-scalar four-fermi coupling A& and the strong cou¬ 
pling a Sj fc, where k sets the momentum scale. Both are 
rising rapidly towards the infrared for momentum scales 
k < 4 — 10 GeV, for see Fig. 13 In this regime chi¬ 
ral symmetry breaking and confinement is triggered and 


takes place in QCD at vanishing magnetic field. Switch¬ 
ing on the magnetic field increases the relevant momen¬ 
tum scale k 2 oc eB and hence decreases A and a s . The 
condensate still grows with B as the R-enhancement in 
the broken phase is still present, only T c decreases. 

Our results from the analytic approach to the quark 
gap equation, presented in Fig. [9] support these findings. 
The position of the turning point R mi „ in both the full 
numerical as well as the analytic analysis of the gap equa¬ 
tion depends crucially on the magnetic field and temper¬ 
ature dependence of the quark gluon vertex, see Fig. [9] 
When contrasted with the quantitative FRG results of 
a s in m, the strong coupling in ( |A2| decays consider¬ 
ably slower towards the UV. In turn, the couplings in the 
qualitative FRG study for full QCD, [74; have a steeper 
decay, for the four-fermi coupling see Fig. |13| Seemingly, 
this already explains the large value of R rniI1 in the cur¬ 
rent DSE-study as well as the small value of R min in [43] , 
which uses approximations similar to m- Note however, 
that we have used the quenched quantitative a s in the 
analytic DSE-study which agrees well with the numerical 
DSE result for \ B ss 1. 

In summary we have identified the physics mechanisms 
behind the T — B phase diagram from our full QCD cal¬ 
culations. Moreover, Fig. [14| suggests a turning point for 
eB min « 1.5— 10 GeV 2 , the large regime for eB min being 
related to the exponential dependence on the couplings. 
Evidently, the effects observed depend on a sensitive bal¬ 
ance of different scales and parameters. Hence, further 
studies are required to fully uncover the intricate underly¬ 
ing dynamics. Very recent findings in AdS/QCD models, 
ES], indicate an inverse magnetic catalysis behavior up 
to eB ss 4 GeV 2 , which supports our findings. 


IV. CONCLUSIONS 

We have investigated the chiral phase structure of 
QCD at finite temperature in the presence of an external 
magnetic field. Our study resolves the discrepancy be¬ 
tween recent lattice and continuum calculations at mag¬ 
netic fields below 1 GeV 2 , see also [43] . We confirm the 
inverse magnetic catalysis effect seen in lattice studies at 
small B. At larger B we see that magnetic catalysis is 
restored, with T c oc V eB. Indications for the turnover 
behavior have already been found in [43], and in [52] 
within two-color lattice-QCD. We hope that further lat¬ 
tice calculations in full QCD at the scales discussed here 
will become feasible soon. 

The reason for this non-monotonous behavior are 
screening effects of the gauge sector, i.e. modifications 
of the gluon self energy, as well as the strong coupling 
a s in the presence of magnetic fields. Moreover we have 
investigated the nature of the chiral transition at finite 
magnetic field. 

Apart from the R-dependence of the critical tempera¬ 
ture, we observe that the phase transition in the chiral 
limit turns smoothly into a crossover with rising B. No- 
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tably, we find a non-degeneracy in the phase transition 
which is due to the explicit isospin breaking caused by 
the different electric charges of up and down quark. This 
non-degeneracy might lead to phenomenological conse¬ 
quences in experimental studies of the QCD phase di¬ 
agram with non-central heavy-ion collisions, as there 
might be a mixed phase between the up and down quark 
transitions. Recent lattice calculations [59j support the 
possibility of a non-degenerate chiral phase transition. 

In addition, our calculations show that, due to this 
isospin breaking, there is a step-like behavior in the up 
quark condensate triggered by the chiral transition of the 
down quark. While this is an significant effect in the chi¬ 
ral limit it smoothens out rapidly with increasing cur¬ 
rent quark mass. Physical current quark masses are in 
the transition regime, and this effect might have phe¬ 
nomenological consequences. To our knowledge, this is a 
novel effect in the QCD phase diagram and it certainly 
deserves further investigation. 

We have used analytic studies of the quark gap equa¬ 
tion and the dynamics of the four-fermi coupling for an 
investigation of the physics mechanisms behind (inverse) 
magnetic catalysis. The results are discussed at length in 
the previous Section |III C[ leading to a rough prediction 
of the turning point at e H min ss 1.5—10 GeV. Our investi¬ 
gations highlight the rich phenomenology of QCD matter 
in external magnetic fields, which motivates further stud¬ 
ies, e.g. at finite chemical potential, towards more realis¬ 
tic descriptions of matter under extreme conditions. Re¬ 
cent studies m have suggested even richer QCD phase 
structures in the presence of magnetic fields. 
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Appendix A: Gluon Propagator and Quark Gluon 
vertex from Dyson Schwinger studies 


Here we discuss the truncation scheme for the quark 
gap equation and the gluon propagator, based on ED ESI. 
The quark gluon vertex is taken as = Zqgqj 11 , with 

WQ 2 ) =d 2 + Q2 ( A1 ) 


Q 2 / A)Ct(/r) log Q 2 /A 2 + 1 

+ A 2 + Q 2 \ 4t r 


25 

5 

(A2) 


containing the parameters 

ch = 7.9 GeV 2 d 2 = 0.5 GeV 2 , 

6 = -18/88, A =1.4 GeV. (A3) 

Here the scales must be identified correctly in order to 
capture the correct dependence with T and eB. We take 
Q to be the symmetric momentum Q 2 = (q 2 + p 2 + 
( q — p) 2 )/ 3 at the vertex with Q 2 = Q 2 + Qq + Q±, 
where Qq = (2ttT) 2 if Qq < (2nT ) 2 and Q\ = 2|eB| if 
Qj_ < 2 1 eB |. We note that this roughly corresponds to an 
identification of scales as in section |III A| with Ab ~ 1, al¬ 
though the present vertex is clearly more sophisticated as 
it includes momentum dependencies and thereby generic 
eB effects. For a current overview of the quark gluon 
vertex in Dyson-Schwinger truncations see IHUSg. Fur¬ 
thermore in order to be able to solve the gluon Dyson- 
Schwinger equation we rely on lattice input for the Yang- 
Mills part, which we then "dress" with magnetic field 
effects, as described above. The reliability of this trunca¬ 
tion was already discussed in detail at finite temperature 
ED and utilized in the presence of magnetic fields before 
|Q. The lattice fit is given by 


ah? 

Q P f Poa(p) logQ 2 /A 2 + 1 
A 2 \ 47 t 

with 

A = 1.4 GeV, c= 11.5 GeV 2 , 

Pa = HV c /3, 7 = -13/22, (A5) 

where a(p) = 0.3 and a and b are temperature dependent 
parameters, which can be found in ES- As discussed 
before the Dyson-Schwinger truncation scheme can be 
related to the skeleton expansion done in our analytic 
estimate, which was motivated by renormalisation group 
invariance 

4 ro ,( 0 »)n,(^)^ B zYM ^ +n/ ,„, (A6) 

where the sum over different polarisation tensor compo¬ 
nents is implied. The right hand side actually serves as 
the input to our numerical study, while the different com¬ 
ponents of n are determined dynamically from solving 
the gluon Dyson-Schwinger equation. 



z YM m — 


Q 2 A 2 


(Q 2 + A 2 ) 2 L VQ 2 + 


Appendix B: Magnetic field dependence of the 
four-fermi coupling from QCD 

As we have discussed in Section IIIIBI the value of the 
NJL coupling A at the intrinsic cutoff scale of the model 
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is determined by QCD dynamics. At large scales the dy¬ 
namics of A is driven by the rightmost diagram shown 
in Fig. El Within simplifications we will motivate the 
functional dependence of this diagram on temperature 
and the magnetic field. In the lowest Landau level ap¬ 


proximation the quarks are constraint to the t-z plane 
denoted by (||), whereas the gluons propagate in all four 
dimensions (||,_L). We write the gluon box diagram in 
Fig. [TT] at zero external momentum as 


f c 

F as {eB> 0.3GeV) ~ 4.5el? 

Jo 


dq\\ , 


«ii 


% 


■ m 2 + a s eBc q 


dq ±, 


q .l 


[q± +qf\+m 2 A + eBa s c A } 2 


(Bl) 


r 


where a s is given as Eq. (12). 


For eB < 0.3 (Bl) is smoothly (quadratic fit) extrap¬ 


olated to e B = 0 with minimising the e !?-dependence. 
The flavor, color and Dirac tensor indices have been con¬ 
tracted, and the comparison with the results for A in 
quenched QCD shown in Fig. [Till si lows that the prefactor 
resulting from the tensor contract is approximately 4.5. 
We have written the propagators in a semi-perturbative 
form with medium dependent mass terms. Further we 
have taken m A ~ 1 GeV as the decoupling scale, m q ss 
300 MeV in the chiral broken phase and c A = c q = 1. 
Strictly speaking both masses are larger than 1 GeV as 
we have to add the cutoff masses oc A 2 . We have chosen 
smaller masses in order to also potentially have access to 


the infrared domain k 0, where the constituent quark 
mass is of the order 0.3 GeV and the gluonic mass gap is of 
the order 1 GeV. Furthermore we have approximated the 
Matsubara sum by an integration, due to the small level 
spacing compared to the magnetic field. This approxi¬ 
mation does not hold small eB, but (Bl I is only used for 


eB > 0.3 GeV. Eq. (Bl I includes the correct dependence 


on a s as well and thus captures eB and T effects qual¬ 
itatively. The model parameters in Section |IIIB| allow 
us to reproduce the quantitative behavior of the chiral 
transition temperature and a more elaborate version of 
Eq. (Bl l does not give much greater insight. Apart from 
the agreement with the T c results from lattice calcula¬ 


tion, Fig. 13 shows that quantitatively reliable results 
from QCD-flows in quenched QCD m are reproduced. 
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